Two-electron quantum dot molecule: 
Composite particles and the spin phase diagram 
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We study a two-electron quantum dot molecule in a magnetic field by the direct diagonalization 
of the Hamiltonian matrix. The ground states of the molecule with the total spin 5* = and 5=1 
provide a possible realization for a qubit of a quantum computer. Switching between the states is 
best achieved by changing the magnetic field. Based on an analysis of the wave function, we show 
that the system consists of composite particles formed by an electron and flux quanta attached to 
it. This picture can also be used to explain the spin phase diagram. 

PACS numbers: 73.21. La, 71.10.-w 



The progress in semiconductor technology has opened 
a rich field of studies focused on the fundamental 
electron-electron interactions and quantum effects in arti- 
ficial atoms and molecules 0]. The most striking feature 
of two-dimensional semiconductor quantum dots (QD) 
and quantum molecules (QDM) is that the correlation 
and magnetic field effects are greatly enhanced compared 
with their normal counterparts. This results in a rich va- 
riety of phenomena in lateral QDMs that have recently 
been investigated experimentally and theoretically, see 
for example Refs. |2[ ||, ^, |^. Also the system param- 
eters can easily be changed, unlike in real atoms and 
molecules where the parameters are natural constants. 
The controllable parameters make it possible to tailor 
the semiconductor structures and for example to switch 
between different ground states. 

In this Letter we concentrate on a two-electron QDM 
consisting of two laterally coupled QDs. In addition to 
the interesting and fundamental correlation and quantum 
effects, this system is very important as a candidate for 
the gate of a quantum computer A central idea is 
to use the total spin {S) of the two-electron QDM as a 
qubit. The ground state spin of the QDM can be either 
5 = or 5 = 1. One of the aims of this Letter is to study 
the regions of 5 = and S = 1 states as a function of 
the important system parameters in their realistic range, 
beyond the approximations of Ref. ^ where a change in 
S was found when the magnetic field B was varied. We 
find a similar crossing at weak B, but also a reappearance 
of the S — Q ground state at larger B. This is caused 
by the correlation effects that are treated accurately in 
our many-body approach. An interesting finding is that 
the system naturally consists of composite particles of 
electrons and flux quanta as in the composite-fermion 
theory 1^. 

We model the two-electron QDM by a 2D Hamiltonian 
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use the one of Ref. 0, namely ^m*ujQ min{(a: — d/2)'^ + 
y'^, {x + d/2)^ -I- y^}. This potential separates to two 
QDs at large inter-dot distances d, and with d it 
simplifies to one parabolic QD. We use GaAs material 
parameters m* /rrie — 0.067 and e — 12.4, and the con- 
finement strength HijJo — 3.0 meV. A is the vector poten- 
tial of the magnetic field (along the z axis) taken in the 
symmetric gauge. One should note that the Hamiltonian 
is spin-free, and the Zeeman-coupling Ez = g*iiBBSz 
(with g* = -0.44 for GaAs) of the magnetic field to 
can be taken into account afterwards. The eigenstates of 
the single-particle part of Eq. (^ are easily obtained as 
expansions 
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where and Uy are integers, and ai is a complex coef- 
ficient. We have used the unit of length as in Ref. [1). 
Fig. |l| shows examples of the non-interacting charge den- 
sities, and Fig. displays single-particle energies. It 
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FIG. 1: Confinement potential Vc (dotted line) and non- 
interacting single-particle density \tjj\'^ along y = for d = 
40 nm. The potential is in units of hujo = 3 meV, and the 
maximum value of each is scaled to unity. The confine- 
ment potential is parabolic in the j/-direction. One can see a 
localizing effect of large d and strong B. 



where Vc is the external confining potential, for which we 



is interesting to compare the localizing effect of B with 
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FIG. 2: Six lowest non-interacting single-particle energies as 
a function of B. The dotted lines are the normal Fock-Darwin 
states {d — 0), and the solid ones are for d = 20 nm. One 
can see that for d > 0, the energies are shifted down and that 
the level crossings occur at weaker B. At some of the level- 
crossings a gap opens, suchs as the one for B « 1.75 T, E m 
7.4 meV. This is due to the lower symmetry of the problem 
for d > 0. In general, there is a reasonable similarity between 
the two sets of energy levels. For sufficiently large d, the 
energies would converge to degenerate states of isolated dots. 
All states are on the lowest Landau level after B ~ 3 T. 

the experimental findings of Brodsky et al. . They see 
a clear splitting of the QDM electron droplet to smaller 
fragments by a strong B. As they work in a low-density 
limit (weak confinement), impurities can have a similar 
role as the potential minima in Fig. |^. The localization 
is also related to the formation of Wigner molecules in 
QDs H], which happens in the low-density limit. 

Similarly to the single-particle states, the full many- 
body wave function with total spin S can be expanded 
as 

^'s(ri,r2) = ^ajj {0,(ri)(/)j(r2) 

+ (-l)^^.(r2)</>,(ri)} (3) 

which is symmetric for S* — and anti-symmetric for 
5=1. Notice that the spin part of the wave func- 
tion is not explicitly written, and we work with spin- 
independent wave functions also below. The coefficient 
vector oti and the corresponding energy Ei for the lih 
eigenstate are found from a generalized eigenvalue prob- 
lem where the Hamiltonian and overlap matrix elements 
can be calculated analytically. Details of the computa- 
tional procedure will be published elsewhere As the 
basis functions we have used all the states with both rix 
and Uy < n with n = 6, and we have checked the conver- 
gence by varying n for all values of d. 

The energy difference AE between the lowest S = 
and S — 1 states is plotted in Figs. and ^ The con- 
vergence of AE can be seen in Fig. lirfor large d. One 
can see that for weak magnetic field values, the ground 
state has 5 = 0. We first concentrate on this regime. For 
B = the 5 = state remains lower with arbitrary d, 




FIG. 3: Energy difference between singlet and triplet states 
AE as a function of the magnetic field B and inter-dot sepa- 
ration d. The white lines separate the 5* = {AE > 0) and 
5 = 1 ground states. 
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FIG. 4: AE for a fixed d = 26.736 nm The effect of 

the Zeeman energy Ez is also shown. The difference between 
expansions with n = 5 and 6 is only around 1 /leV. One can 
see that the second S = state disappears for even a weak 
Ez. 

see e.g. Ref. [|j for a discussion of this exact property. 
However, there is a strong decrease in AE as a function 
of d. This can be understood from the fact that very dis- 
tant QDs interact only weakly and the energy difference 
of the two spin states is smaller. There is also a strong 
decrease in AE as a function of B in the 5 = state. For 
d = 0, AE is rather linear until the crossing to the 5=1 
state, but for d > the curves are rounded. The cross- 
ing point of the different ground states does not depend 
strongly on d; it changes only from 1.6 T to 1.2 T as one 
moves from d = to 40 nm. Due to this, changing the 
total spin in an experimental setup is not easy by just 
changing d. On the other hand, around AE = 0, the 
slope of AE is rather large and changing 5 by S is the 
most natural choice. One can achieve a change of 5 also 
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by changing the strength of Vc- This changes the ratio 
between the energies resulting from the confinement and 
electron-electron interaction. For a weak Vc, the interac- 
tions are stronger and the transition occurs at weaker B 
value. Thus the change of the Vc can be seen as a change 
of the effective value of B. 

The transition from the weak-_B S — state to the 
5 = 1 state is most simply explained by the fact that the 
energies of the two lowest single-particle states approach 
each other as B is made stronger (see Fig. |^). At some 
point this difference is smaller than the exchange-energy, 
and the system spin-polarizes. One can see that also for 
this state, AE decreases strongly as a function of d. 

There exists a second region of 5 = ground state 
around B ^ 6 T at d — 0. The question whether 
this state terminates at large d remains open, as the 
energy differences at d > 40 nm are smaller than the 
error made in the expansion. The existence of this 
S = state for small value of d can be understood 
on the basis of a parabolic two-electron QD, which can 
be shown to have the exact wave function of the form 
^' = {xi2 + iyi2)'"/(f'i2)e~2(''i+''2)^ where m is the rela- 
tive angular momentum and / is a correlation factor [ pl . 
The simple form is due to the separation of the center- 
of-mass and relative motion of the electrons. For B = 
the ground state has m = (and 5 = 0), and when B is 
made stronger, the ground state m has increasing posi- 
tive integer values. These transitions happen because the 
larger m states have smaller Coulomb repulsion between 
the electrons, and as the single-particle energies group 
together to form the lowest Landau level, it is favorable 
to move to larger m. The second S = region in Fig. ^ 
corresponds to m = 2. It is surprising that the second 
5 = region extends to such a large d. For example, at 
B = 4.5 T and d = 35 nm the two QDs of a QDM are 
rather decoupled (see Fig. ||), but still the 5 = state 
remains the ground state. To analyze the structure of 
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FIG. 5: |*[(a;,i/), (d/2,0)]p for the state S = at B = 4.5 T 
and d = 35 nm. The contour spacing is uniform. We mark 
with a plus the fixed electron on the right-hand QD. One can 
see that the electron density is localized to the left QD. There 
is only a small deformation from circular symmetry. Notice 
that the effective inter-dot distance is much larger than d due 
to the Coulomb repulsion of the electrons. 

this state for small and large values of d, we have located 



the vortices of ^. This can be done by finding the zeros 
of \1/ and studying the change in the phase of \1/ in going 
around each of the zeros. A surprising finding is that for 
both large d and d = there are two vortices at both elec- 
tron locations, see Fig. || for an example. The particles 
are thus composite particles of an electron and two flux 
quanta, in similar fashion as in the composite fermion 
(CF) theory The most remarkable feature of the 

finding is the stability of composite particles against the 
change of d. 
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FIG. 6: |*[(x,jy),(d/2,0)]p for the state S = 0atB = 6T 
and d = 10 nm. The arrows (rotated for clarity) depict the 
two flux quanta at the electron fixed at a; > minimum of V^. 

We have done a similar analysis for the 5=1 state 
at B = 3 T and various values of d. We found that the 
many-body state again consists of composite particles, 
but this time there is only one flux quantum per electron. 
One should note that for an odd number of flux quanta 
per electron, the wave function changes sign when the 
particles are exchanged, corresponding to 5 = 1. Simi- 
larly, for even flux numbers the state has 5 = 0. 

If one expands jvE'p for small ri2, one obtains 

1*1' « + -^rfT^' + 0{rl^+^) , (4) 
m+ 2 

where C is the scaled strength of the Coulomb interac- 
tion [H , and m > is the number of flux quanta per elec- 
tron. One can see that for larger m values, the density 
grows more slowly as a function of ri2 (see Fig. ^). One 
should note that the same expansion is valid also for the 
larger electron numbers, resulting from the cusp condi- 
tion 1^ . Thus the electron-electron interaction is smaller 
when the number of flux quanta per electron grows. 

One can use the CF-type approach to explain the phase 
diagram of Fig. When one moves from _B = to 
stronger values of B, the number of flux quanta in the 
system increases. At the first 5 = state there is no 
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fluxes in the system, and at the transition point the flux 
number changes to one per electron. For weakly coupled 
QDs this transition happens at smaller B than in the 
strongly coupled ones, because a distant zero of the wave 
function increases the kinetic energy less than a close one 
does. In the following transition points the number of 
flux quanta changes again by one per electron. The rea- 
soning for the d-dependence of the first transition applies 
to other ones also, and this can be seen in Fig. ||. 

An interesting prospect resulting from the discussion 
above is to use the CF-type approach to describe the 
many-body states of electrons at strong B in various con- 
fining potentials. One should note that after the phase 
structure of the wave function is fixed, one is left with the 
bosonic part of the wave function. The quantum Monte 
Carlo techniques are especially useful for obtaining this 
part 

If one adds the Zeeman term to the Hamiltonian, the 
energy of the S — 1 state is lowered by ~ 25 fieV/T and 
the 5 = energy is unaltered, see Fig. ^. This makes 
the high-B 5 = state to terminate at d value around 
5 nm. One should note that it is possible to lower the 
Zeeman term also in the experimental setup by applying 
a tilted magnetic field. Most probably the singlet-triplet 
separation in energy for i? « 6 T is too small for a qubit, 
but in principle it is possible to change 5 by varying the 
Zeeman term in this parameter range. 

In conclusion, we have determined the total-spin phase 
diagram of the two-electron quantum dot molecule as a 
function of the magnetic field and the inter-dot distance. 
Our results support the possibility to use the system for 
a gate of a quantum computer H. In addition, we have 



found that the system consists of composite particles of 
electrons and the attached magnetic field fiux quanta. 
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